Shear viscosity and spin diffusion coefficient of a two-dimensional Fermi gas 
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Using kinetic theory, we calculate the shear viscosity and the spin diffusion coefficient as well as the 
associated relaxation times for a two-component Fermi gas in two dimensions, as a function of tem- 
perature, coupling strength, polarization, and mass ratio of the two components. It is demonstrated 
that the minimum value of the viscosity decreases with the mass ratio, since Fermi blocking becomes 
less efficient. We furthermore analyze recent experimental results for the quadrupole mode of a 2D 
gas in terms of viscous damping obtaining a qualitative agreement using no fitting parameters. 



INTRODUCTION 



The properties of 2-dimensional (2D) Fermi systems 
are fundamental for our understanding of a wide range 
of phenomena including organic and high-T c supercon- 
ductors, 2D nano-structures, and 3 He films. A new gen- 
eration of experiments are now probing the many-body 
properties of atomic Fermi gases in 2D traps [IH3] . This 
provides a unique possibility to systematically explore 
the physics of 2D systems using the high experimen- 
tal control characterizing atomic gases. Recently, there 
has been a lot of interest in the transport properties of 
atomic gases. One reason is that transport coefficients 
provide excellent probes for strong correlations, since 
they can change by orders of magnitude due to interac- 
tions. It has been shown experimentally that 3D atomic 
gases may form a perfect fluid with a shear viscosity r\ 
having the least possible value consistent with quantum 
mechanics [3]. This has inspired a lot work investigat- 
ing the connections between the physics of atomic gases, 
and other strong coupling systems including quark-gluon 
plasmas, and liquid Helium [5]. Also, recent experiments 
demonstrate that the spin diffusion coefficient approaches 
a scale set by quantum mechanics for a 3D resonantly in- 
teracting atomic gas [6]. The first experiments probing 
the collective mode spectrum of a strongly interacting 2D 
Fermi gas were recently reported [J]. The frequency of 
the breathing mode was shown to be provide evidence of 
a classical dynamical scaling symmetry [8], whereas the 
damping of the quadrupole mode was used as a measure 
for the shear viscosity of a 2D Fermi gas. 

We calculate the shear viscosity and the spin diffusion 
coefficient as well as the associated relaxation times for 
a two-component Fermi gas in 2D using kinetic theory. 
The dependence of the viscosity on the mass ratio of the 
two components is analyzed, and we show that the mini- 
mum value is reduced for systems with a mass imbalance. 
We furthermore analyze the recent experimental results 
for the quadrupole mode in terms of viscous damping [7 
obtaining a qualitative agreement. However, our analy- 
sis shows that further work is needed to understand the 
experiments quantitatively. 



FORMALISM 

Consider a 2D gas of two fermionic species a = 1,2 
with mass ra a and density n a = fep CT /47T so that the 
total density is n — n\ + n-i. The range of the interac- 
tion is taken to be much shorter than the interparticle 
spacing, and there is therefore no interaction between 
identical fermions. We shall focus on two steady state 
non-equilibrium situations: one with a spatially varying 
local mean velocity u(r), and one with a spatially vary- 
ing magnetization M(r) = n\(t) — 112(1"), which we for 
concreteness take to have the forms u(r) = [u x (y), 0] and 
M(r) = M(x). As a result of the velocity field u x (y), 
there is a net current H xy in the y-direction of momen- 
tum along the x-direction, and likewise M(x) induces a 
net magnetization current jm in the ir-direction. Within 
linear response, we can write 

= ~r]dyU x and j M = -Dd x M (1) 

which defines the shear viscosity 77 and the spin diffusion 
coefficient D. 

We briefly outline a variational method to calculate the 
shear viscosity and the spin diffusion coefficient within 
kinetic theory. Further details are given in Refs. [9]- 
IT2] . In kinetic theory, both coefficients are obtained 
from a steady state solution to the Boltzmann equa- 
tion. In the hydrodynamic limit, the distribution func- 
tions /<r(r, p) are close to the local equilibrium form = 
l/[exp(/3^°) + l] with&° = p 2 /2m <y -u(r)-p-p <7 for the 
case of a local velocity field u and = p 2 /2m a — (r) 
appropriate for a local magnetization M. Here /ii(r) and 
fi2 (r) are the spatially varying chemical potentials corre- 
sponding to the magnetization and ft = 1/T (we use 
units where fcs = h = 1). When these local equilibrium 
functions are plugged into the left side of the linearized 
Boltzmann equation, it can be written as 

9 fn d ^ = I« and d -f^=h (2) 
oe ay Oe Ox 

where l a is the collision operator for component a, and 
f„ = l/[exp/?(e CT — Ha) + 1] is the equilibrium function 
with e a — p 2 /2m a . Here $JJ = p x p y /m a and $^ = 
p x jm a for shear viscosity and spin diffusion respectively. 
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The momentum and spin currents are given by 

n xy = J d 2 k{^h + $2/ 2 ) 

3m 



J d 2 fc($f f x $f h) 



(3) 



with d 2 k = d 2 k / (2ir) 2 . To proceed, we need an approxi- 
mate solution to the Boltzmann equation In 3D, the 
ansatz 8f a oc ^/"(l — /°) for the deviation of f a (r,p) 
away from equilibrium is known to yield results within 
2% of the exact result for the viscosity [H [13] . We there- 
fore use this ansatz for the 2D case which yields 



77 = 2/3- 



($»i 2 ) 2 



and D 



(<S> d2 ) 2 



X (<S> D H[<i> D } 



(4) 



as variational expressions for the viscosity and spin dif- 
fusion coefficient. We have defined the average (<I> 2 ) = 
2- 1 Eal^kfM-fM and* = d{ ni -n 2 )/d^-^) 
is the magnetic susceptibility. The linearized collision in- 
tegral can after symmetrization be written as 



d 2 k x d 2 k 2 



Pr 



dO 



, da 
~dW 



x(A^/i U / 2 u (l-/i u )(l-/2°) 



(5) 



where p r = (maPi — ra\p 2 )/M is the relative momen- 
tum of the incoming scattering particles, 9' is the angle 
between the outgoing and incoming relative momenta, 
and da/dO is the differential cross section. The to- 



- m2, and m~ = + 
We have defined the function 
$1 + $2 — $3 — $2 f° r shear viscosity and 



$3" + for spin diffusion. It de- 



tal mass is M = mi 
is the reduced mass. 
A$ 

A$ = $f 

termines the contribution of a given collision to the mo- 
mentum and spin transport respectively. The reason for 
the factor 2 difference in the expressions for r\ and D in 
Q is that the two components contribute with the same 
sign to the momentum current and with opposite signs to 
the magnetic current, see This also causes the sign 
differences in the expressions for A$ in the two cases. 



Relaxation times 

The viscous and spin relaxation times r n and Tp> which 
give the typical time between collisions for the two types 
of motion are useful for estimating whether a system is 
in the hydrodynamic regime. Suitable definitions can be 
obtained by writing the collision integral as la- — 5f a /r. 
which gives r) = 2T n f3((^) 2 } and D = t d /3({<S> d ) 2 )/x- 
Performing the integrals yields 



V 
n 



fo°°deef 



and 



D 
n 



TP 

2mx 



(6) 



where we have taken mi = m 2 and n± = n 2 = n/2 
for simplicity. This reduces to rj = nT n ep/2 and D = 
cftd/™- in the degenerate limit, whereas 77 = ut^T and 
D = TpT/m in the classical limit. We have used x — 
to/2tt for T <C Tp and x — n/2T for T > Tp. 



Scattering cross section 

When the range of the interaction is much shorter than 
the typical interparticle spacing, the scattering between 
the (7 = 1 and a = 2 fermions is predominantly s-wave. 
The 2D cross section for relative momentum p r is a = 
m 2 \T (p 2 / 2m r )\ 2 / p r with the T-matrix given by [T4HT6] 



T(e) 



2tt 



1 



m r ln(|£^(,|/e) + in 



(7) 



which has a pole at a 2-body bound state with energy 
E b = —l/2m r a 2 . 

CLASSICAL LIMIT 

Consider the classical limit T 3> Tp a = k F(J /2m a . In 
this limit /" « 1, and the integrals in Q are straight- 
forward to perform. We obtain for the viscosity 



T 



2n 2 



n\n 2 I V (T/E b ) 



(8) 



with 



2 



[ln(|£ 6 |/T)-0.92] 2 + 7r 2 



(9) 



The shear viscosity depends only on the reduced mass 
in the classical limit. This is because the scattering 
only depends on the relative coordinates in this limit, 
since there is no Fermi-blocking. For fixed total density 
rii+rt 2 j the viscosity is minimum for n\ = n 2 as expected, 
since the scattering becomes less frequent with increas- 
ing population imbalance. The viscosity in the classical 
limit for mi = m 2 and n± — n 2 was reported while this 
manuscript was being written |17j and the result agrees 
with (|8|-([9]) for that case. 

Likewise, Q yields for the spin diffusion coefficient in 
the classical limit 



T 



An 2 nI D {T/E b ) 



(10) 



with 



I D (T/E b ) = / die 
Jo 



ln(\E b \/Tt) 2 +7r 2 



1 



[ln(|£ 6 |/T)-0.42] 2 +7r 2 



(11) 
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Here, we have for simplicity taken wi\ = m-i and n\ = 
n,2 = n/2 and used \ — njlkoT in the classical limit. 

Equations ([8])-([9]) and (10|-([TT| should be compared 
with the analogous expressions obtained for the 3D 
case: r\ = l5{mk B T) 3 / 2 /32^ [TTJ CEB] and D = 
3 v ^(fc B T) 3 / 2 /16 v ^ : [12 for a classical gas in the uni- 
tarity regime. The reason for the more complicated In- 
dependence in 2D is the intrinsic energy dependence of 
the T matrix (JT]), which means one never recovers the 
simple power law predictions for an energy independent 
cross section: r\ ~ nvml m f cx \T and D ~ vl m f oc yT 
where l m { ~ 1/na is the mean free path. 



NUMERICAL RESULTS 

In this section we present numerical results for the vis- 
cosity and spin diffusion coefficient obtained from Q. 



Viscosity 

In Fig. [l] (a) , we plot the viscosity as a function of T 
for the mass ratios mi/7712 = 1 and m\jrai = 6/40. The 
latter corresponds to a mixture of 40 K and 6 Li atoms. 
We have taken the density of the two components to be 
equal and \n(\Eb\/Tp 2 ) = 1- For high temperatures, the 
viscosity approaches the classical result (J8|, whereas it 
increases strongly for low T due to Fermi blocking |19j . 
This results in a minimum of the viscosity at T ~ 0.6Tp2 
for ?7ii/m2 = 1, whereas the minimum is located at 
T ~ 1.3Tp2 for mi/m2 = 6/40 due to the larger Fermi 
temperature for the light component a = 1. 

An important result is that the minimum viscosity 
of the mass imbalanced mixture is significantly smaller 
than for the mass balanced mixture; for the mass ra- 
tio TO1/TO2 = 6/40 it is a factor 0.6 smaller. This can 
be understood as follows: changing both masses keep- 
ing mi/rri2 = 1 clearly does not reduce the minimum 
value of rj/n, since this simply amounts to rescaling Tp; 
however, reducing mi while keeping mi fixed makes the 
Fermi blocking less efficient on the scale of Tp\ and the 
minimum value of rj/n is reduced essentially since the 
classical result ^ holds for lower T/Tpx- The mini- 
mum value of the viscosity is subject to intense inter- 
est due to a conjecture inspired by results for a certain 
class of strong coupling theories ;20j, which states that 
the ratio of the viscosity over the entropy of any sys- 
tem obeys the universal bound -q/s > 1/Att [3T]. In the 
inset of Fig. [I] (a), we therefore plot 77/s for the same 
parameters as in the main plot. The entropy density 
s = s\ + s 2 is obtained from the ideal gas expression 
s a = -Jd 2 k[f°lnf° - (1 - /°)ln(l - /»)]. Again, we 
see that the minimum value of rj/s is significantly smaller 
for the mass ratio m\jm,2 = 6/40. Intriguingly, it seems 
to follow from kinetic theory that a two-component sys- 




5.5 
5 



CM 

I 35 

^ 3 



1.5 







(b)/ 










rr^/m^eMO 











| 1 og(IE h l/ EF J 



FIG. 1: (color on-line) (a) The viscosity as a function of tem- 
perature for \n(\Eb\/Tw 2) = 1 and mx/rtii = 1 (blue line), 
m\/m2 = 6/40 (cyan dash-dot line). The classical limits 
are plotted as red dashed and purple dotted lines. The inset 
shows 77/s as a function of T. (b) The viscosity as a function 
of \n{\Eb\/TF2) for T = Tf2 and mi/ma = 1 (blue line), and 
mi/m2 = 6/40 (cyan dash-dot line). 



terns with a sufficiently large mass ratio can break the 
conjectured bound rj/s > I /Air. A similar effect is in fact 
present for 3D systems. 

In Fig.[T](b), we plot rj/n as a function of In(|i?h|/Tp2) 
for T = Tp2- The viscosity is minimum in the strong 
coupling regime \n(\Eb\/TF 2 ) ~ 0{1) as expected. In the 
classical limit, it follows from ^ that the minimum is lo- 
cated at \n(\E b \/T) ~ 0.92. With decreasing m 1 /m 2 , the 
minimum moves to larger values of hi(\Ei > \/Tp2) because 
Tpi increases. 
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FIG. 2: (color on-line) The spin diffusion coefficient as a func- 
tion of temperature for ln(\Eb\/TF2) = 1 and mi/mz = 1. 
The classical limit is plotted as a red dashed line. The inset 
shows D as a function of the interaction strength for T — Tf-2 ■ 



Spin diffusion coefficient 

The spin diffusion coefficient is plotted in Fig. ~~| 
as a function of temperature for m\jm 2 = 1 and 
\n(\Eb\/Tp 2 ) = 1. For high T, it approaches the clas- 
sical value ( 10 ) whereas Fermi blocking makes it in- 
crease strongly for low T, leading to a minimum value 
at T = 0.85Tp. The inset shows D as a function of 
ln(\Ei,\/Tp 2 ) for T = Tp2- Again, the minimum value 
is for \n{\E b \/T F2 ) 



0(1). For high T, ((TO) predicts 
0.42 to be the minimum value. 



Validity of kinetic theory 

Let us briefly discuss the range of validity of kinetic 
theory. For weak coupling | ln(|~~&|/e.F|)| 3> 1, the kinetic 
approach is accurate except for extremely low tempera- 
ture, as has been shown for the viscosity in 3D [18]. For 
strong coupling, one must expect the Boltzmann equa- 
tion ~~|) to break down for low temperature where there 
are no well-defined quasi-particles. In 3D, calculations of 
the viscosity based on the Kubo-formalism show that the 
kinetic approach is accurate down to temperatures signif- 
icantly below Tp, even for strong coupling p~3j . n~8l . l22j . r23] . 
We expect a similar result to hold in 2D. In particular, 
kinetic theory is likely to be reliable at the tempera- 
tures where we predict 77 and D to be minimum. The 
occupation of the closed channel molecule can further- 
more have significant effects on thermodynamic proper- 
ties in 2D [~~~ . Similarly, corrections to the single channel 
approximation for the T-matrix ~f|) could influence the 
transport properties considered here, although one would 



expect small effects for a broad resonance. 



EXPERIMENTS 

The frequency and damping of the quadrupole mode 
of a 2D Fermi gas of 40 K atoms with equal populations 
in two hyperfine states were recently measured ~fj. The 
results were interpreted in terms of viscous damping ap- 
propriate for the hydrodynamic regime. The amplitude 
damping of a collective mode can be calculated from ~~5] 



|(£mech)t| 
2(-E'mech)t 



(12) 



where (-Emech)t is the time averaged mechanical energy 
of the mode. Taking the velocity field of the quadrupole 
mode to have the form u(r) = (x, —y) cosut, we get 



(^moch)t = y / d 2 rn(r)u 2 {r), 



(13) 



where we have used that the potential energy of the mode 
is equal to the kinetic energy, and we have neglected any 
interaction energy. The viscous damping is for this ve- 
locity field, following ~~5]~~~], given by 



(E 



mcch, 



= -2b 2 



d 2 r 



(14) 



Here we have used the real part of the complex dynamical 
viscosity t](cj) = r]/[l — iuiT v (r)] ~~7] evaluated at the 
quadrupole frequency ujq to obtain a cut-off in the outer 
classical regions of the cloud, where the viscosity is given 
by ~~|) with mi = m 2 and n\ = n 2 = n/2, and therefore 
is independent of density. In the classical regime, ( 12 ) 
becomes using © and (6) 



r d = 



2uj 



irVNI v 



du- 



(15) 



where N — N\ + N 2 is the total number of particles 
trapped and to is the 2D trapping frequency. 

In Fig. [3] (a) , we plot the damping of the quadrupole 
mode taking T/T F = 0.47 and N = (Ep/huj^) 2 = 4300 
particles trapped which are the experimental parameters 
appropriate for Fig. 1 in Ref. "~J. We have calculated 
the damping as a function of \n(kpa 2 ) and the x's are 
the experimental results reported in Ref. ~~j. We see 
that the theory agrees qualitatively with the data. In 
Fig. ~5] (b), we plot the damping as a function of T for 
various coupling strengths taking N — 3500 particles 
trapped to model the experimental situation of Fig. 3 
in Ref. ~~j. Again, the theory accounts qualitatively for 
the experimental results which are plotted as x 's. Note 
that we have no fitting parameters. For both sets of data, 
the agreement between theory and experiment is best for 
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FIG. 3: (color on-line) (a) The damping of the quadrupole 
mode as a function of interaction strength. The x 's are 
the experimental results of Ref. [7|. (b) The damping of 
the quadrupole mode as a function of T for various coupling 
strengths with the x's the experimental results of Ref. [7J. 
The inset shows the viscous collision rate oj±t v . 



large \n{kpa2) and large T/Tp whereas there are signif- 
icant quantitative discrepancies in the strong coupling 
regime of small \n(kpa2). Similar results were reported 
in Ref. [T7] for the spatial average of the viscosity using 
the classical limit approximation. 

It is perhaps surprising that the agreement is best in 



the weak coupling regime where the system is collisionless 
rather than hydrodynamic. This is illustrated in the inset 
in Fig. 3 (b) which shows uj^t^: the hydrodynamic condi- 
tion Lu±T ri < 1 is fulfilled only for \n(kpa,2) — 2.7, where 
as cu±T n > 1 for larger ha.(kpa2) indicating collisionless 
dynamics. However, despite being based on hydrody- 
namics, the viscous damping approach turns out to work 
rather well in the collisionless regime. In 3D it has in fact 
been shown to yield exact results in the collisionless limit 
provided one uses the complex dynamical viscosity eval- 
uated at the collisionless collective mode frequency EH). 

The reason for the discrepancy between theory and 
experiment for small ln(fc^a 2 ) where the system is hy- 
drodynamic, can be strong coupling effects making the 
kinetic approach quantitatively inaccurate as discussed 
above. Better agreement could also be obtained by solv- 
ing the Boltzmann equation approximately by taking mo- 
ments with basis functions for 6f a ESESEH]- Such an 
approach has indeed been successful in describing the fre- 
quency and damping of collective modes in 3D. 

To summarize, we have, using kinetic theory, calcu- 
lated the shear viscosity and the spin diffusion coeffi- 
cient for a two-component Fermi gas in 2D. Both trans- 
port coefficients have a minimum value somewhat below 
the Fermi temperature. We showed that the minimum 
value of the viscosity can be reduced significantly with 
increasing mass ratio of the two components. Using a 
viscous damping approach, we qualitatively accounted 
for recent experimental results for the damping of the 
quadrupole mode of a 2D Fermi gas. However, our re- 
sults also showed that further work is needed to obtain a 
quantitative understanding of the results. 

I am grateful to M. Kohl for discussing the experimen- 
tal results in Ref. [7] with me. 



[3] 
[4] 
[5] 

[6] 



[7] 

[8] 
[9] 

[10] 

[11] 
[12] 



B. Frohlich et al. 



Prog. Phys. 72, 



[1] K. Martiyanov, V. Makhalov, and A. Turlapov, Phys. 

Rev. Lett. 105, 030404 (2010); 
[2] M. Feld et al, Nature 480, 75 (2011): 
Phys. Rev. Lett. 106, 105301 (2011). 
A.Sommer et al., |arXiv:1110.3058| 
C. Cao et al, Science 331, 28 (2011). 
T. Schafer and D. Teaney, Rep. 
126001(2009). 

A. Sommer et al., Nature 472, 201 (2011); A. Sommer, 
M. Ku, and M. W. Zwierlein, New J. Phys. 13, 055009 
(2011). 

E. Vogt, M. F eld, B. Frohlich, D. Pertot, M. Koshorreck, 
and M. Kohl |arXiv:1111.1173| 
J. Hofmann, [a"rXiv:1112.1384p 

H. Smith and H. H0jgaard Jensen, Transport Phenomena 
(Oxford University Press, 1989). 

G. Baym and C. J. Pethick, Landau Fermi-Liquid Theory 
(Wiley, 1991). 

P. Massignan, G. M. Bruun, and H. Smith, Phys. Rev. 
A 71, 033607 (2005). 

G. M. Bruun, New. J. of Phys. 13, 035005 (2011). 



6 



[13] G. M. Bruun and H. Smith, Phys. Rev. A 75, 043612 
(2007). 

[14] S. K. Adhikari, Am. J. of Phys 54, 362 (1986). 
[15] D. S. Petrov and G. V. Shlyapnikov, Phys. Rev. A 64, 
012706 (2001). 

[16] M. Randeria, J.-M. Duan, and L.-Y. Shieh, Phys. Rev. 

Lett. 62, 981 (1989). 

[17] T. Schafer, |arXiv:1111.7242"l 

[18] G. M. Bruun and H. Smith, Phys. Rev. A 72, 043605 
(2005). 

[19] D. S. Novikov, arXiv:0603184. 

[20] G. Policastro, D. T. Son, and A. O. Starinets, Phys. Rev. 

Lett. 87, 081601 (2001). 
[21] P. Kovtun, D. T. Son, and A. O. Starinets, Phys. Rev. 

Lett. 94, 111601 (2005). 



[22] T. Enss, R. Haussmann, and W. Zwerger, Ann. Phys. 

326 770 (2011). 
[23] S. Riedl et al, Phys. Rev. A 78, 053609 (2008). 
[24] J. P. Kestner and L.-M. Duan, Phys. Rev. A 74, 053606 

(2006) ; W. Zhang, G.-D. Lin, and L.-M. Duan, ibid 77, 
063613 (2008). 

[25] L. D. Landau and E. M. Lifshitz, Fluid Mechanics (Perg- 

amon Press, 1987). 
[26] G. M. Bruun and H. Smith, Phys. Rev. A 76, 045602 

(2007) . 

[27] T. Nikuni and A. Griffin, Phys. Rev. A 69, 023604 (2004). 
[28] S. Chiacchiera, T. Lepers, D. Davesne, and M. Urban, 
Phys. Rev. A 84, 043634 (2011). 



